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We study the symmetry properties of the Hubbard model with spin-orbit interactions of Rashba 
and Dresselhaus type. These interactions break the rotational symmetry in spin space, so that 
the magnetic order cannot be excluded by using the Bogoliubov inequality method. Nevertheless, 
we rigorously show that the existence of the magnetic long-range orders may be ruled out when 
the Rashba and Dresselhaus coupling constants are equal in modulus, whereas the 7/-pairing can 
be always ruled out, regardless of the microscopic parameters of the model. These results are 
obtained by imposing locally the SU(2) gauge symmetry on the lattice, and rewriting the spin-orbit 
interactions in such a way that they are included in the path ordered of the gauge field on lattice. 
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I. INTRODUCTION 

One of the most interesting features of the interacting 
systems is the existence of a macroscopic order which may 
break the underlying symmetry of the Hamiltonian. For 
instance, it is known that the continuous rotational sym- 
metry in the three-dimensional spin space of the isotropic 
Heisenberg ferromagnet is broken by the spontaneous 
magnetization that exists in the limit of vanishing mag- 
netic field for a three-dimensional latticed For systems of 
restricted dimensionality and at any finite temperature, 
it has been argued that there is no macroscopic order, 
both on the basis of heuristic arguments and by a rigorous 
approach based on the use of an operator inequality due 
to Bogoliubov,- The piece of work by Bogoliubov stim- 
ulated numerous investigations on the upper and lower 
bounds for thermodynamic averages, stating that sponta- 
neous ordering is usually not possible in low-dimensional 
systems £ 

When a lattice chain of spins is considered, pioneering 
work by Mermin and Wagner- ruled out magnetic order- 
ing both in one-dimensional (ID) and two-dimensional 
(2D) quantum spin systems with an isotropic (Heisen- 
berg) interaction between the spins, at any nonzero 
temperature. There are, however, other interactions, 
such as the Ruderman-Kittel-Kasuya-Yosida interaction 
(RKKY) j£ in which the conduction electrons mediate the 
interaction between spins. These interactions are long- 
ranged and hence not covered by the Mermin- Wagner 
result. Theoretical work has shown that ordering is not 
possible in ID RKKY systems^ but nothing similar has 
yet been attempted for the 2D case. Nevertheless, re- 
cently it has been shown that the Mermin- Wagner theo- 
rem can be applied also to systems of lattice spins which 
are spin coupled to itinerant and interacting charge carri- 
ers^ Within this approach, it has been rigorously proved 
that neither (anti-)ferromagnetic nor helical long-range 



order is possible in ID and 2D, at any finite tempera- 
ture. This conclusion is not related to a specific model 
Hamiltonian but applies to a wide class of models in- 
cluding any form of spin-independent electron-electron 
and single-electron interactions. 

Nevertheless, this argument becomes inconclusive 
when interactions which explicitly break the spin sym- 
metry are considered. In this case magnetic order cannot 
be excluded. This is the case for instance of spin-orbit 
interactions (SOIs) such as the Rashba^ or the Dressel- 
haus^ ones. We recall that the Rashba spin-orbit field 
can be attributed to an electric field that originates from 
a structural inversion asymmetry whereas the Dressel- 
haus interaction comes from bulk inversion asymmetry. 10 
Since both these spin-orbit interactions break the sym- 
metry of the spin interaction, their presence means that 
the magnetic order cannot be excluded. For complete- 
ness we mention that this conclusion is valid except for 
the special case in which the Rashba and Dresselhaus 
contributions are of equal magnitude,- This result turns 
out to be experimentally relevant because the Rashba 
spin-orbit interaction can be electrically tuned to match 
the Dresselhaus interaction, opening up possibilities for 
tuning magnetic order electrically. 

Here, we investigate a specific model, namely the 
single-band Hubbard model, in presence of SOIs, and 
we rigorously show that the existence of long-range or- 
ders may be ruled out also when SOIs are present and 
the Bogoliubov inequality method is applied. Since SOIs 
break the spin symmetry of the model, we first introduce 
SOIs in the model by imposing local SU(2) gauge sym- 
metry on the lattice, then we rewrite SOIs in such a way 
that they are included in the path ordered of the gauge 
field on the lattice. In this way, the SU(2) invariancc 
is restored allowing the use of the Bogoliubov inequality 
and the derivation of upper bounds for order parameters. 
At this stage, it is worth mentioning that local symme- 
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tries, such as the local gauge symmetry, cannot be bro- 
ken spontaneously. Indeed, the Elitzur's theorem-^ states 
that a local gauge symmetry cannot be broken sponta- 
neously, that is, the expectation value of any gauge non- 
invariant local observable (order parameter) must vanish. 
This result means that the spontaneous breaking of the 
gauge invariance can only occur when the local symme- 
try is explicitly broken by gauge hxing. Hence, first one 
chooses a gauge and in this gauge the remaining global 
gauge symmetry is spontaneously broken, as it happens 
for instance in the Anderson- Higgs mechanism,— We also 
mention that for two-dimensional systems, if the residual 
symmetry is continuous, after the gauge fixing, at finite 
temperatures the associated order parameter must vanish 
according to the Mermin- Wagner theorem. Nevertheless, 
when a discrete symmetry is considered then the ordered 
phase can still exist. These results will be used in the 
next sections to describe the Hubbard model with SOIs 
within a SU(2) lattice gauge theory. 

We point out that, recently the model here investi- 
gated has been extensively studied. Indeed, the analysis 
of the properties of non-centrosymmetric superconductor 
CePtsSi stems from the single-orbital Hubbard model 
with SOIs, since this model has been considered as a 
valid effective model for describing low-energy quasipar- 
ticles in the Fermi liquid stated. Moreover, experiments 
of the penetration depth of I^PdaB and I^PtaB reflect 
the difference of the strengths of the SOIs. Motivated by 
these data, the study of pairing symmetry of these non- 
centrosymmetric superconductors is again based on the 
single-band Hubbard model with SOIsJ^ New develop- 
ment of technology of epitaxial growth makes it possible 
to fabricate heterointerface between different transition- 
metal oxides. In SrTiOa/LaAlOs the mobile electrons 
at the heterointerface are mainly introduced to 3d or- 
bitals of Ti 3+ ions. Considering the crystalline electric 
field by a ionic model, the two-dimensional d xy orbital 
has the lowest orbital energy level in 3d manifold, as 
also shown by x-ray absorption spectroscopy with the 
linearly polarized light. To microscopically model these 
hybrid structures, a two-dimensional single-band Hub- 
bard model with the SOIs as a minimal model of the 
two-dimensional electron gas at the heterointerface has 
been considered^ Finally, we mention that a purely the- 
oretical study of the interplay between electron-electron 
interactions has been recently performed again within a 
single-band Hubbard model with on-site, and eventually 
nearest-neighbor matrix elements of the Coulomb inter- 
action in the presence of SOIs. 16 

The paper is organized as follows: in the next section 
we summarize some results on the lattice gauge theory, 
emphasizing their application to the Hubbard model pre- 
sented in Sect. III. Section IV describes the application 
of the results of the previous sections and the main result 
of the paper, namely the application of the Bogoliubov 
inequality to the Hubbard model in the presence of SOIs. 
The last section is finally devoted to remarks and con- 
clusions. 



II. REVIEW ON HAMILTONIAN LATTICE 
GAUGE THEORY 

To straightforwardly understand the outcomes of the 
next section, here we summarize some important results 
of the lattice gauge theoryJi Since the interactions may 
be interpreted as a consequence of the local invariance 
of the theory under some gauge transformation groups, 
our aim is to show that the SOI may be included in the 
Hubbard model as the path ordered of gauge fields on 
the lattice. 

If G is a gauge transformation group, its elements are 

g = e ie aT% 

where 0\,...,6 n are the group parameters and T\,...,T n 
are the group generators. These operators fulfill the con- 
dition 

[T ai Tb] — ifabcTc, 

where (fabc)2bc=i are the group structure constants. 

In quantum mechanics, the group transformations are 
realized by means of unitary operators F in the Hilbert 
space V s 

| a >->-| a' >= F \a>, 

for any \a >£ V s . Therefore, if A and ^(x) are the 
crystal lattice and a field operator, and x € A, then 

*(ar) F + V(x)F = e i6aTai S{x), 

where the operator F is 

In this expression Q a — J2 x eA^ + ( x )'^a^{x) with a — 
l,...,n are the so-called charge operators and they are 
a representation of the group generators, thus imply- 
ing that they verify the same commutation relation. If 
[Q a ,H] — 0, then the quantum system has the global 
symmetry under the group considered. 

A local transformation is defined by the following 
relation: 

F+y(x)F = e Wa ^ T "*(x). (1) 



From these preliminary definitions, it is obvious that a 
system is globally gauge invariant if the model Hamilto- 
nian contains only terms like ^> + (x)^(y). On the other 
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hand, to get the local gauge invariance of the model we 
introduce a link gauge field U(x, x + e;), where e; is a 
unit vector along the I direction, with I = 1, ...,d, and a 
link vector potential A a (x) = A a (x,x + e/). Therefore, 
for each pair of (x,x + e/) nearest-neighbor lattice sites 
we write: 

*(x)+*(x + ej) * + (x)J7 ; (x)*(x + ej), 

with 

C/i(x) = C/(x,x + e ; ) = e^^ )r -. 

If lattice sites (x, y) are not nearest-neighbor sites, we 
consider the link path ordered sequence 

n-1 

r(xi,x„) = (J (xjjXi+i), 

1=1 

and we define U(x, y) as 

U(x,y)= Yl U(z,z + ei) = e ij: ^-'^^ A '^ z)Ta , 

where the link gauge field is such that 

F+U(x, y)F = e ieJx)T -U(x iy )e^ e ^ y)T -. (2) 



Now, we introduce the link group generators E l a (x) as 
follows 

d 

F = JJY[ e i8 a (x)(*+{x)T a H!(x)+E l a (x)) = Y[ e iOa(x)G a (x) ^ 
x£Al=l i£A 

(3) 

where 

[G a (x),G b (x')]=if abc G c (x)5 x , x ,. (4) 



Therefore, if [G a {x),H] = then the Hamiltonian H 
exhibits the local gauge symmetry. 

As a special case, if F is an infinitesimal transforma- 
tion, by means of Eqs. ([I])-© we get the commutation 
rules 



[G ,*(x)] = -T Q *(x), (5) 
[G a ,*+(x)} = *+(x)T a , 



[G a ,Ui(x)] = -[T a , Ui(x)], 

with 

G a = ^ G a (x), 

x£A 

and for each pair of lattice sites (x, y) 

[G a {x),U(x,y)]^T a U{x,y) 1 (6) 
[G a (y),U(x,y)] = -U(x,y)T a . 



III. HUBBARD MODEL WITH SOI 

After having considered the main results for the lattice 
gauge theory in the preceding section, let us now discuss 
how they can be applied to the Hubbard model in the 
presence of SOIs. 

It is well-known that the Hubbard model describes 
the electronic motion in the crystals within the tight- 
binding approximation. 18 If one introduces the spin 



operator \&(x) = 




, where c CT is the annihila- 



tion operator at site x 6 A c 1 d for an electron with 
spin cr, then the Hubbard Hamiltonian may be written as 

H = E t(x-y)iS! + (x)*(y) + Uj2Mx)n(x)> (?) 

x : y£A x£A 

where n a (x) = c+(x)c CT (x) is the particle number 
operator with spin a. This Hamiltonian is rotationally 
invariant in the spin space, i. e. [H, S a ] = 0, where 
S a = ExeAS a (x) with 

S + (x) = c^(x)c l (x) S~(x) = c|(x)ct(x), 

S z (x) = ±(n t (x)-n l (x)). (8) 

The operators defined by Eq. ([SJ are the SU(2) group 
generators, and their commutation rules are: 

[S + ,S~] = 2iS z [S*,^] = ±S ± . (9) 

According to Mermin- Wagner theorem, a spontaneous 
magnetic order is absent in the Hubbard model^ This 
result has been obtained adding to this Hamiltonian 
of SU(2) symmetry-breaking term, and applying the 
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Bogoliubov inequality. We notice that the SU(2) sym- 
metry exhibited by the Hubbard model is a global one. 
To get a SU(2) local gauge symmetry for this model, 
we may implement the lattice gauge theory previously 
introduced. In this way the Hubbard Hamiltonian can 
be written as 



H= J2 t(x-y)^ + (x)U(x,y)^(y)+Uj2 r H( x ) n i(^)- 



with an appropriate Lande factor g^, and the spatial 
components are defined as 



A 3 {x)=0, 



A\x) 



2mc I P \ X2i \ 2mc / o 
a , Af{x) = — I P 



eh 







eh 







For instance, in the continuum the substitution 



(14) 



(10) 



where U(x, y) is a lattice gauge field and H gaU ge is the 
free gauge field Hamiltonian. 

Introducing the following field operators 



p 1 e ^a a , 2 . a0 a 



-> — (p- -A a — y 

2m 2m c 2 



2 " 



allows to obtain the well-known Hamiltonian with SOI 
Rashba and Dresselhaus couplings in the presence of a 
magnetic field: 



r = J a (*) = £(£ Ela (x) + s°0«0), 



x£A 



i. e. the transformation generators of the local SU(2) 
group, the Hamiltonian Eq. (|T0|) is invariant under the 
considered group, that is [J a {x),H] = 0. This result can 
be easily verified using Eqs. ©-([H]), that for SU(2) group 
are 



Here. 



and 



H=^+H R + H D+ g L ^.B(x). 



H R = -(P x a 2 - Pya 1 ) 



(15) 



[J a ,U(x,y)] =-[—,U{x,y)l 



[J a (x),U(x,y)]=—U(x,y), 



(11) 



H D = l j[P x o l - P y a 2 ). 

are the Rashba and Dresselhaus terms, respectively. 

Now, let us introduce SOI in Eq. (J?} by using the 
SU(2) lattice gauge theory previously outlined. To this 
end, we define for our case the link gauge fields as 



[J a (y),U(x,y)} 



~U(x,y) Y , 



and 

[J a Mx)} 



(12) 



[J a ,*+ (x )]^+ (x )^_, (13) 



where a 1 are the Pauli matrices. 

It is well-know that SOIs in condensed matter physics 
can be described in terms of a non-Abelian gauge the- 
ory^ by means of a non-Abelian potential A^(x) — 
%Afa) (a = 1,2,3) where 



eh 



B a (x) A ia {x) 



eh 



e^E {x) 



Ui(x) = e^ i ^ aaAla ^ = e i 9iPcr 1 +^ 2 ) 



A i I ■ i\ 2 (cos(f><T 1 + i singer 2 ) 



Ai A a e-*^-T) 



U 2 (x) = e~ i ^ aaA2a(x) = e - i 9{^ 1 +P<y 2 ) 



Ail — i\ 2 (singer 1 + i cosher 2 ) 



Here, B l {x) are the components of the external magnetic 
field and E l (x) is the electric field produced, for example, 
by nuclei in molecules or solids, and a, i, j = 1,2, 3. 

For 2D models with SOIs interaction of Rashba and/or 
Dresselhaus type, the time component of SU(2) potential 
is 

A a0 (x) = ^^ B B a (x), 



Ai -A a e** 
A 2 e"^ Ai 

U 3 (x) = 1, 



with g 



(16) 



Ai = cos(gV a 2 + f3 2 ), A 2 = sm(gy/ a 2 + f3 2 ), 
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and tan <j) = % . 

We notice that the link gauge fields verify the following 
commutation rules: 

[A 1 (x),A 2 (x')] =2^(a 2 -/3 2 ) f 7^ 

[Ui(x), U 2 {x')] = 2i(A 2 ) 2 (cos 2 <b - sin 2 4>)<r 3 . (17) 

It is worth stressing that when the coupling constants a 
and (3 are equal, the symmetry group becomes an abelian 
C/(l) group. Furthermore, within the definition reported 
in Eq. (fTl)) . the gauge fields U[(x) are independent on 
the choice of the lattice site. This means that, if the 
pair of lattice sites {x, y) in the hopping term of the 
Hamiltonian Eq. ([TU|) are not nearest-neighbor sites, 
then U (x, y) depends on sites x and y only. But this is 
not the general case and U(x, y) may depend on the link 
sequence connecting the x site to the y site. Therefore, 
if r = (x, xi, x n —i, y) is a path on the lattice, then 
U r (x,y) = U il (x)Ui 2 (x 1 )...U in {x n ^i) where i k is the 
direction on the lattice connecting x^-i and x k of the 
path r. We point out that it is possible to obtain the 
general form of U {x, y) by imposing the time reversal 
symmetry of the Hamiltonian in Eq. (1101) . The time 
reversal operator is given by 

xeA 

where K is the complex conjugation. By imposing 
T+HT = H one obtains 

U(x,y)=( 9f?'V\ J { X / y \). (18) 
v ' y > \-j*(x,y) g*{x,y) J v ' 

Here, the functions g(x,y) and j(x,y) are connected to 
the electronic jump from x to y with spin conservation 
and spin flip, respectively. These functions depend obvi- 
ously on the electronic path followed in the jump. 
Assuming the local SU(2) gauge symmetry on the lat- 
tice, the interaction is contained in the path ordered of 
the gauge field on the lattice U(x,y), so that the group 
SU(2) algebra can be used when we apply the Bogoliubov 
inequality to get the upper bound on the order parame- 
ter. 



IV. APPLICATION AND RESULTS 

It is important to observe that the spin-orbit coupling 
breaks explicitly the SU(2) spin symmetry, and there- 
fore we cannot say anything on the spontaneous mag- 
netization. Nevertheless, we will show that within our 
approach, if Rashba and Dresselhaus coupling constants 
are equal, then the Hubbard Hamiltonian in Eq. (|10p ex- 
hibits a C/(l) rotational symmetry in spin space and the 



SU(2) can be restored by a gauge transformation in the 
spin space. This conclusion will be used to prove that if 
hopping matrix t(x — y) ^ 0( r~^ra ) then the magnetic 
ordering is absent in d = 2, for any finite temperature, 
in agreement with Mermin- Wagner theorem. In d = 1, 
it is always possible to restore the SU(2) symmetry in 
spin space by a gauge transformation, and therefore the 
magnetic order is equally absent. On the other hand, the 
r] pairing superconductivity is vanishing for any value as- 
sumed by a and /?, because the U(l) symmetry in pseu- 
dospin space is not broken at all. To better clarify these 
points we consider separately the study of the magnetic 
order in the d = 2 and the d = 1 cases, and then the r\ 
pairing long-range. 

To rigorously look at the absence, for finite tempera- 
tures, of the spontaneous long-range orders mentioned 
above, one should first remove the degeneracy of the 
model Hamiltonian and then study the expectation val- 
ues involved. The degeneracy is usually removed by in- 
troducing a symmetry-breaking term to the Hamiltonian 
under study. Then, the main steps to follow are: a) 
the identification of the order parameter involved in the 
phase transition; b) the addition to the Hamiltonian of 
a symmetry-breaking term relevant for the order under 
study; c) the adequate choice of the operators in the Bo- 
goliubov inequality to single-out the order parameter; d) 
the search for non trivial upper bounds; e) the proof that, 
in the thermodynamic limit and in the low-dimensional 
cases, the order parameter vanishes, at any non-zero tem- 
perature. The procedure now outlined will be applied in 
the next subsections to rule out the possibility of mag- 
netic long-range order and then to the superconducting 
case. 

A. Magnetic order in two dimensional lattice with 

a = ±[3 

Here, we show that in the special case when Rashba 
and Dresselhaus SOI become equal in intensity, the mag- 
netic order is excluded in two-dimensional lattices. When 
Rashba and Dresselhaus coupling constants are such that 
a — ±/3, the link gauge fields, Eqs. (IT5|) . are: 

Ui(x) = U X = Ail + i^(±a 1 + a 2 ), 
v2 

U 2 (x) = U 2 = Ail -z^=(a 1 ±a 2 ). (19) 
v2 

They are commuting operators \U\, U 2 ] = 0, as it can be 
deduced looking at Eq. (|17p. Thus, the symmetry group 
is the U(l) abelian group of rotations around a given 
direction. This statement implies that the ordering of 
the link gauge fields is not important and we can write 
the hopping term in the Hamiltonian in Eq. (Jjj) as follows 

t(x-y)V + (x)U(x,y)y(y) = 

x,y£A 
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= t ( x - {x)U 1 {x)U 1 {x + ei)* 

:r,y(z A 

*...*Ui(x+lei)U2(x+lei)U2(x+lei+e2)...Ui(x+lei+me2)* 
**(i + Ie 1 + me 2 ), (20) 

where y = x + le,i + me2, &i and e 2 being the unit vectors 
on the X\ and X2 axis of the two dimensional lattice, 
respectively and 1 and m are integers. Moreover, we can 
also write 

U(x,y) = Uf^Uf-* 2 = e is'('-™)«±-* (21) 

with g' = g\J a 2 + ft 2 and n± = -^(1,±1,0), where the 
sign + (-) holds if a = j3 (a = —(3). It is easy to verify 
that if a = /3 then the Hamiltonian in Eq. ([7]) is rotation- 
ally invariant around the ri+ = -i (1,1,0) axis, that is 

[n- S,H] =0 where S is the total spin. This is also easily 
observable in the Hamiltonian (|15[) . whereas if a = — /3 
the symmetry is around the ri_ = i(l, — 1,0) axis. Af- 
ter straightforward algebra, the gauge interaction can be 
written as 

U(x,y) = e -VMz)«±^ e «?'My)™±-^ (22) 

where 8± (x) — \f 7 ln± ■ x. 

Defining a gauge transformation F as 

p _ J~J e i9's T (x)n ± -s(x) ^ 

we find that the Hubbard Hamiltonian with spin orbit 
coupling is globally SU(2) invariant in the spin space: 

F + V + (x)U(x, y)V(y)F = *+(x)*(y). (23) 

If we consider the average value of the spin projection 
along £3 axis, we can write 

< S 3 (x) >= tr(e-P H S 3 (x)) = tr(F+e' 0H FF+S 3 (x)F) = 

= tr(e-? F+HF p-S(x)), 

where pi = Ri^e^ji and R is a rotation matrix of 0^(x) 
around n± axis. The Hamiltonian is now SU(2) in- 
variant, and then by the global gauge transformation 

O = e -ia'e^{x)a ± -s we get 

< S* 3 (:z;) >=< S 3 {x) >', 

where < >' is the average value of the S 3 (x) in the 
transformed Hamiltonian. 

Since for the Hubbard model without spin orbit cou- 
pling the spontaneous magnetization is absent^ intro- 
ducing the magnetic order parameter as 

m h (Q, A) = i- ]T e iQx < S 3 (x) >= JL < S 3 (Q) >', 
' ' xeA ' ' 

we trivially deduce 

lim lim mniQ, |A|) = 0. 

h->0 |A|-K» 



B. Magnetic order in one dimensional lattice 

We consider the Hubbard model on a chain of lat- 
tice sites with a constant lattice a assumed for simplicity 
equal to 1. Therefore, the crystal lattice is A C Z, and 
the Hubbard Hamiltonian is 

|A| |A|-i |A| 

fl = EE *0')* + (0*(i+i) + tf£ n W»U0- 

(=1 j=l 1=1 

In one spatial dimension, we have one link direction, so 
that the lattice may be seen as a sequence of links. Also, 
it can be considered as U(l) subgroup of SU(2) where we 
have one gauge field associated with the link (1,1 + 1), 
that is, the group is abelian and its elements are: 

Ui(x) = e~ i ~^ a "' A}a = Ail + iA 2 (cos(/)cr 1 + singer 2 ). 

To introduce the spin-orbit coupling we perform the re- 
placement 

# + (0*(/+j) -> * + (l)U 1 (l)U 1 (l+l)...U 1 (l+3-l)y(l+j) = 

= *+(i)tf(M + j)*(/ + j). 

Since the gauge fields are independent on the site, we 
write 

U(l, l+j) = e ia ' jH - s , (24) 

with g' = g\/u 2 + /3 2 and n = (cos <f>, sin</>, 0). 

Summarizing, the Hubbard Hamiltonian is rotation- 
ally invariant around the vector n. Indeed if S = 
Si='i "5(0 denotes the total spin operator, then we have 
[S- n, H] = 0. Furthermore, Eq. (|2~4")) is very similar to 
Eq. (|2T|) . so that we may define the unitary operator 
F = ni=i e i9 ' ln - S M. Therefore, we can write 

F + V+(l)U{l, I + + j)F = *+(/)*(/ + j), 

that is, the transformed Hamiltonian is the Hubbard 
model without the spin-obit interaction. Therefore, using 
the argument previously outlined for the two-dimensional 
lattice, we can conclude that the spontaneous magneti- 
zation is absent for any a and /3 values, also in the ID 
case. 



C. 77 pairing long-range order in one and two 
dimensional lattices 

The 77 pairing long-range order is connected to the 
spontaneous symmetry breaking of the £7(1) invariance 
around the third axis in the pseudospin space of the 
Hamiltonian. 22 The Hubbard Hamiltonian exhibits this 
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symmetry and the introduction of the spin orbit inter- 
action does not modify this property. If we define the 
pseudospin operators as 

rf{x) = Cf(x)ci(x) 77+0) = c^(x)c+(x), 



r] (x) = -(n-fO) + H(x) - 1), 



then [if,H] = 0. Here, 



Then, from the commutator 

[B,*+(x)U(x,y)*(y)} = (e^ x -e^y)^+(x)U(x,y)^(y), 
we get the average value of the double commutator as 
< [B+,[B,H]] >= 

= 2 ^2 t(x - y)(l - cos(q- (x - y)))* 

x.y£A 



From the commutators 

[r] 3 ,^{x)] = [?7 3 ,* + (x)] = (25) 

it is easy to prove that 

[r, 3 ,^+(x)U(x,y)^(y)] = 0. 

The absence of long-range r\ pairing in the Hubbard 
model without spin-orbit coupling has been widely stud- 
ied^ Here, following the same approach, we extend this 
result to the Hubbard Hamiltonian with SOI. To this end, 
we introduce the symmetry breaking external field A in 
the Hamiltonian as follows 

H^H-X( v + (Q)+ v -(-Q)), 

and we define the rj pairing order parameter as 

A(Q) = lim JL £ e -*Q-« < v + {x ) >= 
|AK» |A| ^ 



* < 9+{x)U(x,y)9(y) > -A(< r?+(Q) > + < r?-(-Q) >). 

Now, if one defines the scalar product by (A, B)=< 
A + B >, then by Schwartz inequality one gets 



so that 



< 9+(x)U{x,y)9(y) 2, 



2q 2 



| < [B+, [B, H}] > | < |A|(^- + AA(Q, A)). (27) 

Since < >= |A|A(Q,A) and 

£<L4,A+] + ><|A| 2 , 



by using the Bogoliubov inequality Eq. (|26l) the proof is 
accomplished. Indeed, we get 



|A(Q)| 2 



d d q 



1 



(27r) d q 2 + Af A(Q) ^ pfc s T' 



Thus, by solving the integral we find out that, when 
A 0, A(Q) —> for d=l, 2, at finite temperature. 
This implies that the Hubbard model with SOI does not 
exhibit the 77 pairing long-range order. 



V. CONCLUSIONS 



lim 

|A|->-oo 



<v + (Q) > 

|A| ■ 



To show that long-range r\ pairing is absent in this model 
we have to show that 

lim A A = 0. 
So, in the Bogoliubov inequality 



< [A, B] > | 2 < 



1 



<[A,A+] + ><[B+,[B,H}] >, 



2k B T 

we define the operators A and B as 

A(q)=r)+(q + Q) B{q)=rf{-q). 



(26) 



We presented an extension of the Mermin- Wagner the- 
orem for the Hubbard model in the presence of SOI, and 
showed that spontaneous magnetic order is ruled out in 
two dimensions, at finite temperature, if the Rashba (a) 
and the Dresselhaus (j3) spin-orbit interactions are such 
that a = ±/3. On the contrary, in one-dimension the 
magnetic order can be excluded, regardless of the val- 
ues assumed by the spin-orbit coupling constants. We 
notice that, when Q—0 in m(Q) the ferromagnetic or- 
der is forbidden, while choosing Q in such a way that 
exp(iQi?j) = ±1 when Ri connects sites in the same 
sublattice and different sublattices, respectively, we ar- 
gue that also the antiferromagnetic order is forbidden. 
We also proved the absence of long-range 77 pairing, at 
finite temperatures, in one- and two-dimensions, inde- 
pendently on a and (3 interaction parameters. As stated 
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for the magnetic order, looking at the ^-pairing order pa- 
rameter A(Q), we may infer that for Q = the s-wave 
pairing can be excluded, for Q = ±ir we exclude the 
^-pairing, and finally for Q ^ {0, ±7r} we rule out the 
existence of generalized 77-pairing order with momentum 
Q. 

For copleteness, we note that the Mermin- Wagner the- 
orem follows from the fact that in low-dimensional cases, 
a diverging number of infinitesimally low-lying excita- 
tions is created at any finite temperature, and thus 
the assumption of a non vanishing order parameter is 
not self-consistent. This consideration, as well as the 
rigorous proof, does not apply at T=0, implying that 
the ground-state may be ordered. For instance, two- 
dimensional ferro(anti)magnetism is possible at zero tem- 
perature: quantum fluctuations oppose but do not pre- 
vent the appearance of a two-dimensional magnetically 
ordered phase. In contrast, for one-dimensional systems 



quantum fluctuations become so strong that they usually 
prevent even ground state ordering. Indeed, it is known 
that the ground-state of the one-dimensional Hubbard 
model is a non magnetic singlet at any band filling and 
for any value of Coulomb interaction U. More generally, if 
the energy spectrum has a gap, it can be shown that the 
model under investigation does not exhibit long range 
order, and interestingly, this energy gap plays the role 
of the temperature in conventional Bogoliubov inequal- 
ity3."26 
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